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ABSTRACT 

Old Globular Clusters (GCs) in the Milky Way have ages of about 13 Gyr, placing their formation time in 
the reionization epoch. We propose a novel scenario for the formation of these systems based on the merger of 
two or more atomic cooling halos at high-redshift (z > 6). First generation stars are formed as an intense burst 
in the center of a minihalo that grows above the threshold for hydrogen cooling (halo mass r\j 10^ Mq) 
by undergoing a major merger within its cooling timescale (^150 Myr). Subsequent minor mergers and 
sustained gas infall bring new supply of pristine gas at the halo center, creating conditions that can trigger 
new episodes of star formation. The dark-matter halo around the GC is then stripped during assembly of the 
host galaxy halo. Minihalo merging is efficient only in a short redshift window, set by the ACDM parameters, 
allowing us to make a strong prediction on the age distribution for old GCs. From cosmological simulations 
we derive an average merging redshift (z) = 9 and narrow distribution Az = 2, implying average GC age (tage) = 
13.0 ±0.2 Gyr including ~ 0.2 Gyr of star formation delay. Qualitatively, our scenario reproduces other 
general old GC properties (characteristic masses and number of objects, metallicity versus galactocentric radius 
anticorrelation, radial distribution), but unlike age, these generally depend on details of baryonic physics. In 
addition to improved age measurements, direct validation of the model at z ^ 10 may be within reach of 
ultradeep gravitationally lensed observations with the James Webb Space Telescope. 

Subject headings: galaxies: high-redshift — galaxies: general — globular clusters: general — cosmology: 
theory 


1. INTRODUCTION 

Globular Clusters (GCs) are compact stellar systems with 
characteristic ma ss ^10^ Mq and radius of a few pc (Heg- 
gie & Hut|2003| ). Their ubiquitous presence around galaxies 


and their old stellar populations make them tools to inves¬ 
tigate early star formation, as sembly history of hos t galax¬ 
ies, and co smological models ( |Katz & Ricotti|[2QT3| |Brodie| 
et al.|2014| ). For exa mple, GC ages were used to constrain the 
age o f the Universe ( [Jimenez et al.|1996} [Krauss & Chaboyerj 
2003| ), providing ear ly independent support to the c oncor¬ 
dance ACDM model ( [Planck Collaboration et al.|2013| ). 

Yet, GC formation remains a debated topic with the lack 
of an established scenario matching all observations. Early 
proposals identified t he high-z Jeans mass (^ 10^ Mq) with a 
protoglobular cloud ( [Peebles & Dicke|1968 ). More recently, 
high-z GC formation as been proposed within dark-matter 
(D M) halos of mass Mh ^10^ Mq, with cooling driven by 
H 2 ( [Padoan et al.[| 1997b, as result of shocks induced by a 
reionization front ( |Cen|200T| ), or as purely baryonic systems 
because stream velocity displaced gas from its parent halo 
( [Naoz & Narayan|2014[ ). Other scenarios focus at somewhat 
lower redshift, ranging from cooling-induced fragmentation 


of (proto)-galaxies ( [Fall & Rees||1985|), to formation during 
galaxy merging/interactions ( Ashman & Zept|1992t Muratov[ 


& Gnedi^[2010[ [Li & Gnedin[|2014|) or within hig h-density 


region s of galactic disks ( [Kravtsov & Gnedin|2005[ Kruijssen 

2014[ ). Analogies with t oday’s young massi ve clusters have 
also been proposed (e.g., |Bastian et al.j2Q13[ ). 

The lack of consensus on GC formation may indicate multi¬ 


ple formation mechanisms. In fact, the oldest GCs have essen¬ 
tially uniform ages centered at ^ 12.8 Gyr and spread ^5% 
(comparable to relative age uncertainty), but wide range of 
metallicities: from 1% to 1/3 the solar value with high metal¬ 
licity systems preferentially at small galactocentric radii. In 
contrast, younger systems show a well defined age-metallicity 
anti-c orrelation and higher metallicity at larger galactoce ntric 


radii ( Forbes & Bridges|2010[[Marm-Franch et al.|2009[ ) 

In this Letter we aim at predicting the probability dis¬ 
tribution function of ages for the sub-population of “old” 
GCs starting from DM halo assembly and merging, which 
are fully determined (in the statistical sense) by fixing the 
ACDM parameters. Our approach parallels studies of evo¬ 
lution of galaxy properties with redshift, such as luminos- 
ity/stell ar mass functions and clustering, all linked to halo as¬ 
sembly ([Trenti et al. [20 10b [[Lacey et al.|201 HfThcchella et al. 


[MT^IBehroozi et aL|2013| 


We propose that old GCs formed via major merging of DM 
minihalos {Mdm ^10^ Mq) that are gas rich (no previous 
star formation) and metal-enriched through outflows originat¬ 
ing from nearby (proto)galaxies. First generation stars form 
during the merger-triggered burst, while subsequent gas in¬ 
fall through minor mergers can lead to multiple populations. 
Since the halo merger rate is approximately c onstant per unit 
redshift (both for major and minor events; [Fakhouri et al.| 
2010| ), there is only a short window of opportunity for this 


mechanism. The redshift needs to be low enough (z < 15) so 
that minihalos are relatively common, but not too low (z ^ 6), 
otherwise the merger rate drops. 
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2. GC FORMATION THROUGH MINIHALO-MINIHALO MERGERS 

We propose that old GCs form when a host galaxy like the 
MW is in the earliest stages of its assembly and lacks a well 
defined disk structure. Its Lagrangian region is still a moder¬ 
ate linear overdensity at z ^ 10 containing hundreds of mini- 
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Fig. 1.— Snapshot at z = 10 {tiookback ~ 13.3 Gyr) of the region that by z = 0 
will collapse to the halo of mass 8 x 10^^ Mq shown in the right panels of 
Figure^ A linear overdensity containing hundreds of collapsing/collapsed 
minihalos is present. Regions involved in GC formation through minihalo 
mergers are highlighed in red (past major merger) and yellow (future merger). 
A movie capturing redshift evolution is available in the online edition. 

halos that are being hierarchically assembled, form stars, pol¬ 
lute their neighbours with metals, and merge to build up a 
massive system by z = 0 (Figure [T]). 

We assume that compact star clusters form during major 
mergers of two gas-rich and star-poor minihalos, whose gas 
has been previously enriched by outflows from nearby halos 
toZ>10-2Zo. The combined minihalo mass needs to exceed 
the threshold for HI cooling {T^tr > 10"^ K), but if either of 
the progenitors was above the cooling mass threshold, it must 
not have started forming stars yet. Under these conditions, 
the merger will shock-heat and concentrate the gas, leading 
to efficient formation of a star cluster. Following the merger- 
triggered starburst. Type II supernovae clear the gas remain¬ 
ing in the minihalo (escape velocity is a few tens km/s), sup¬ 
pressing further star formation after ^10 Myr. Subsequently, 
slow winds from AGB stars produce chemically enriched gas 
(jD’Ercole et al.|20081), which is retained in the potential well, 
and diluted by infall of new pristine gas through minor merg¬ 
ers and/or accretion. This gas loading triggers one or more 
bursts of star formation, more centrally concentrated than the 
first generation, with merger mass ratios and timing after AGB 
pollution imprinting diverse chemical signatures in i ndividual 
clusters, consistent with observations ( |Renzini|2008| ). The GC 
minihalo is eventually incorporated into the merger-tree main 
branch, with continued mergers and tidal interactions strip¬ 
ping DM and leaving a “naked” GC by z = 0. Qualitatively, 
the scenario is illustrated in Figure 

3. QUANTITATIVE MODELING 

To evaluate the plausibility of the scenario depicted in Sec¬ 
tion [3 and compute the GC formation rate, their ages, and 
spatial distribution at z = 0 we res ort to a DM-only cosmolog¬ 
ical simulati on run with Gadget2 (|Springel|2005 ) , with setup 
described in |Trenti et al.| ( |2009] ); |Trenti & Shull ( |2010| ), and 


|Ramirez-Ruiz et aL| ( |2Q15| ), but tailored volume/mass reso¬ 
lution and updated cosmology (Ua o = 0.685, Um.o = 0.315, 


UAQ"0»0462,(78 = 0.p 8, Hs = 0.9585, h = 0.673; |PlanckC6r 
laboration et al.|201^ . We simulate a 10^ Mpc^ (comoving) 


volume from z = 150toz = 0 using N = 403^ particles, while 
we run a high-resolution version of the same initial conditions 
to z = 5.5 (tiookback ^12.8 Gyr) with A = 812^ particles (mass 
resolution 8.3 x 10"^ M(^. Our multi-scale approach is based 
on assigning membership of high-resolution particles to z = 0 
halos through mapping of particles IDs, and studying halo as¬ 
sembly at higher redshift from the full resolution run, saving 
snapshots at uniform intervals of 10 Myr. 

The high resolution DM-only simulation is post-processed 
with the star forma tion and chemical enrichment model of 
Trent! et~ar] ( |2Q09| ) to fiag minihalos that are in gas-rich. 


star-free conditions. Specifically, we follow the formation 
of Population-Ill (metal-free) stars in halos above the mini¬ 
mum threshold for H 2 cooling with a time-evolving Lyman- 
Wemer background, triggering Population-Ill star formation 
when a pristine halo crosses the mass M^/„(z) shown in Fig¬ 
ure 1 of [Trent! et al.| ( |2009| ). Under these conditions, a halo 
is flagged as enriched by Population-Ill stars, and then starts 
forming Population II (metal-enriched) stars once above the 
hydrogen cooling limit (T^ir > 10"^ K). As basic description of 
metal pollution within our DM-only framework, Population- 
II halos, once they begin forming stars, are assumed to have 
outfiows propagating with spherical symmetry at fixed speed 
^wind = 60 km/s (appropriate for dwarf-like galaxies). Ha¬ 
los that grow with redshift staying below Mmin(z) until z ^ 13 
(point at which Mminiz) corresponds to atomic cooling) may 
either form a “late-time” Population-Ill star, or a Population 
II cluster/dwarf galaxy if polluted b y outfiows from nearby 
(proto)galaxies dTrenti & Shull|2010| ). 

As zeroth order characterisation of metal enrichment, we 
track for each minihalo the number ^ of nearby halos that have 
metal outfiows active for long enough to reach its center. We 
consider outflow-enriched (pollution counter ^ > 1), star-free 
halos with T^ir ^ lO'^K as potential GC birth-sites. The condi¬ 
tion that needs to be satisfied in our framework to create a GC 
is a major merger with an another smaller outfiow-enrichec0 
minihalo that increases the mass so that T^ir > 10"^ K. If ei¬ 
ther progenitor has T^ir > lO^K pre-merger, then the merger 
must happen before stars are formed. Otherwise supernova 
feedback clears the minihal o of gas. Assuming a typical cool¬ 
ing time of (1 -2) X 10^ yr ( Tegmark et al. 1 1997 ), we define a 
timescale Ar = 150 Myr from crossing of the cooling thresh¬ 
old (T^ir > lO'^K). After such time we no longer consider a 
GC as the outcome of a major minihalo merger. Furthermore, 
we define the minimum mass ratio between progenitor halos 
as M 2 /Ml > 1/4. 

Finally, we save IDs for all particles involved in GC forma¬ 
tion to track halo membership, check for subsequent mergers, 
and characterize their spatial distribution at z = 0 through ID 
re-mapping into the lower resolution run. 

Our simplified gas treatment from DM-only simulations ig¬ 
nores the relative velocity difference between baryons and 
CDM imprinted at recombination time, which can signifi- 


^ This guarantees that GC minihalos are wel l resolved with N > 500 par¬ 
ticles, sufficient to characterise their properties JTrenti et al.|2010a}. 

^ Mergers between a metal-enriched and a metal-free, star-free minihalo 
would also satisfy these requirements but such configurations are rare since 
pairs of merging minihalos have strongly correlated 
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Fig. 2.— Qualitative illustration of GC formation at high-z through a minihalo-minihalo merger (Sec.[^. 


cantly impact formation of high-z halos wit h 10^ K < T^j 


< 


10^ K, suppressing the halo ma ss function ([Tseliakhovich & 


Hiratal 

et al.|| 


l||2010[ |Naoz et al.||2011|), and the gas fraction (|Naoz 
201 3|). Quantifying the impact of stream velocity for 


GC formation is difficult without a hydrodynamic simulation, 
because of competing effects. Gas-rich minihalo mergers re¬ 
quire progenitors where Population-Ill star formation is sup¬ 
pressed. In our framework this happens because of H 2 pho¬ 
todissociation, but stream velo city works as well. In addition , 
the halo bias would be higher ( Tseliakhovich & Hirata|2QlQ| , 
enhancing the merger probability. Overall, the number of gas- 
rich minihalo mergers may be only moderately suppressed, or 
possibly even enhanced, in contrast to the ^ 50% suppression 


of star formation inferred for th e general population of mini¬ 
halos ( |Bovy & Dvorkin|[2013| ). We expect our estimate on 
the GC birth-rate to be accurate within a factor two or better 
compared to a full treatment of stream velocity. 


3.1. Merger rate and age probability distribution 

Figure shows predictions for the GC formation rate in 
the most massive simulated halo {Mk = 4.9 x 10^^ Mq, se¬ 
lected for best statistics). At z > 5.5 (Universe age < 1 Gyr), 
there are 279 minihalo-minihalo mergers fulfilling the con¬ 
ditions for GC formation. The average merger redshift is 
(z) = 9.3 (tiookback = 13.2 Gyr) with dispersion Az = 1.88 (age 
spread ~ 200 Myr). Minihalo-minihalo merging is strongly 
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Fig. 3.— Top panel: Formation rate of GCs through minihalo-minihalo 
mergers in the most massive dark matter halo of our simulation (M/j = 4.9 x 
10^^ Mq at z = 0). Data have been smoothed (boxcar filter width 50 Myr). 
Atz>5.5, A^ = 279 clusters are formed with (z) =9.3 (dispersion Az = 
1.9). Bottom panel: Cumulative number of GCs formed per unit volume 
(Mpc~^ comoving) versus lookback time, for the two most massive halos of 
the simulated volume (solid/dashed blue), for the next three halos in mass 
ranking (solid, dotted and dashed red), and box-averaged value (black). 


peaked with redshift because of competing processes. At 
z > 12 atomic cooling halos are rare, and there is not sufficient 
time for chemical pollution by outfiows originating in nearby 
star forming sit es. Later, as the red shift decreases so does 
the merger rate (|Fakhouri et al.||2010|), and at z < 8 structure 
formation starts shifting toward larger scales. This is quali¬ 
tatively similar to the circumstances leading to the “Gamov 
peak” in Nuclear Astrophysics. 


Figure [^is representative of the average GC formation rate, 
albeit with significant halo-to-halo scatter. Its bottom panel 
illustrates the cumulative GC number for the five most mas¬ 
sive halos in the box comprising a mass range (0.8-4.9) x 
10^^ Mq, as well as the box average. The number of GCs is 
normalised by the comoving volume of each halo in the ini¬ 
tial conditions (proportional to halo mass). The box average 
normalization is the simulation volume. Figurej^depicts halo- 
to-halo variations of a factor two. Furthermore, it shows that 
more massive halos (the five blue and red lines) have slightly 
older GCs compared to less massive halos (the box average 
line) by about 100 Myr. This is a consequence of the ear¬ 
lier chemical enrichment and biased halo mass function that 
is found in overdense environments. 

With an average time delay of ^ 250 Myr, the large ma¬ 
jority of the products of minihalo-minihalo mergers has a mi¬ 
nor merger (progenitor mass ratio > 1/10) or receives a mass 
infall >10^ Mq between two snapshots (10 Myr temporal 
spacing). Major mergers are also frequent: ^ 25% of the first 
generation clusters has one. M ore generally, from Extended 
Press-Schechter modeling (see [Trenti et ar]|2008| ) we derive 
that a halo with Mh = 10^ Mq at z = 10 will evolve into a de¬ 
scendant that 500 Myr later (by z = 5.8) has doubled in mass 
at > 98% confidence. Therefore, the high accretion rate on 
the minihalos suggests that there is sufficient gas to mix, di¬ 
lute chemically the AGB ejecta from first generation stars, and 
concentrate them at the halo center. These con ditions can lead 
to for mation of multiple stellar populati ons ([Carretta et al. 
2009| ), widely observed in galactic GCs (IPiotto et al.||2007| ). 
However, detailed modeling requires hydrodynamic simula- 


tions and is beyond the scope of this initial work, focused on 
making predictions on the robust (parameter-free) probability 
distribution of GC ages. 

Our simulation does not include halos especially represen¬ 
tative of the MW/local group because of its small volume. 
Nevertheless, since minihalo mergers are a universal pro¬ 
cess, we can rescale the results to estimate the number of old 
galactic GCs. The MW halo mass is debated, ranging from 
8 X 10^^ Mq to ^ 4 X 10^^ Mq at 95% (e.g., van der Marel 


et al.|2Q12{|Phelps et al.|2013||Kafie et al.|2014| ), almost en 

compassing the halos shown in Figure!^ In these halos there 
are 279 to 32 old GCs, suggesting fulTconsistency between 
observed frequency and predictions in our framework (old 
GCs in the MW are about 50% of the total population, that 
is A ^ 80). Overall, our theoretical estimate of the old GC 
number is accurate within a factor two, comparable to halo- 
to-halo variations. For improved estimates, GC survival to 
z = 0 is also important, as well as detailed modeling of the hy¬ 
drodynamics (stream velocity, impact of progenitor mass ratio 
and cooling timescale At). 

However, all these baryonic processes, including stream ve¬ 
locity, do not affect to first approximation the predicted age 
distribution, which is the key observable that can falsify our 
model. In this respect, we note that the GC “formation epoch” 
shown Figure is that of the minihalo-minihalo merger, and 
thus an upper limit. Stars form in the following ^ 10^ yr (first 
generation), and possibly a few hundred Myr later when fur¬ 
ther gas accretion triggers subsequent generations. Assuming 
2 X 10^ yr as typical delay for the average stellar age since 
the merger, our best estimate of the absolute age of old GCs 
is ^ 13.0 Gyr with a la spread of ^ 0.2 Gyr. This is in 
full agreement with the current determin ation of ages for the 
old GC population tage = 12.8 ± 0.6 Gyr (Marm-Franch et al. 
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Fig. 4.— Top-left panel: Projected spatial distribution of DM particles in the most massive simulated halo at z = 0 {M^ = 4.9 x 10^^ Mq). Bottom-middle 
panel: Spatial distribution (same redshift and same halo of top-left panel) of the subset of particles that have been part of the N = 279 minihalo-minihalo mergers 
at z > 5.5, which we propose are associated to the oldest GCs. Bottom-left panel: Projection of average value of chemical enrichment tag for the total metallicity 
of GCs (light t o dark for increasing numb er of pollution events), highlighting increasing metallicity at small galactocentric radii in qualitative agreement with 
observations by |Marm-Franch et al.H2009) . Right column contains the same panels but for the fifth most massive simulated halo (M^ ~ 8 x 10^^ Mq). 


2009| ). We predict that, if old GCs formed through mini-halo 
mergers, improvements in absolute and relative age calibra¬ 
tions should converge toward GC formation during the epoch 
of reionization and age scatter as small as ^0.2 Gyr. 


3.2. GC masses, galactocentric distribution and tidal 
stripping 

Our model quantitatively predicts the age distribution of 
GCs formed through minihalo-mergers, an d pr oduces the cor¬ 
rect order of magnitude of objects (Section [3?T] ). To further es¬ 
tablish its plausibility, we investigate the expected GC masses, 
galactocentric distribution and efficiency of tidal stripping in 
removing DM around GCs by z = 0. 

A typical star formation efficiency used in studies of high- 
reds hift objects is =0 .03 (stellar to baryonic mass ra¬ 
tio; I Alvarez et al. 2012| ). This gives for a minihalo with 


Mh ^ 10^ Mq a stellar mass 4 x 10^ Mq, consisten t 

with the average GC mass ^ 10^ Mq (Heggie & Hut 2003), 
even if there is a diffuse stellar component formed during the 
merger (and later stripped). Furthermore, our model has a cut¬ 
off scale given by the HI cooling mass, naturally predicting a 
distribution of masses peaked around a characteristic value 
(log-normal type), as observed for GCs in both the Milky 
Way and exter nal galaxies (e.g. [HarrisjParmentier & 


Gilmore|2005| ). Thus, the model has no need to invo ke pref¬ 


erential disru ption of lower mass star clusters (e.g. [Fall & 
Zhang|20OT]). 


Finally, we discuss model predictions for galactocentric GC 
distribution and efficiency of tidal stripping of the DM enve¬ 
lope. We tag in the z = 0 snapshot the particles associated 
to GC formation in the high-resolution, high-z run. Results 
are shown in Figure tagged particles (middle panel) have 
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Fig. 5.— Density profile versus galactocentric radius at z = 0 for the halo 
shown in the right panels of Figure]^ The DM density (black-solid line) is 
well represented by a NFW profile ^ashed-blue line), while particles asso¬ 
ciated to minihalo-minihalo mergers have a steeper profile (solid red), con¬ 
sistent with p ~ observed in the distribution of Galactic GCs (dotted 
magenta). Profiles have been normalized at r = 10 Kpc h“^. 


higher galactocentric concentration compared to all particles 
(top panel). The z = 0 radial density profile of GC particles is 
shown in Figure for a typical non-interacting halo, demon¬ 
strating that it follows the slope p{r) ^ observed in the 
distribution of MW GCs. This is not surprising sinc e mini 
halo merger s are > 2 a DM peaks at z ^ 9, which Moore 
et al. ( 2Q06| ) found to be distributed with a p{r) ^ profile 


by z = 0. The middle panels of figure also show a diffuse 
distribution of particles involved in mimhalo mergers. This 
supports our proposal that tidal interactions during the merger 
process strip GCs of their initial DM halos, although the re¬ 
sults need to be confirmed by a full high-resolution simulation 
to z = 0. Last, through particle tagging we derive the cor¬ 
relation between z = 0 galactocentric distance and c hemical 
enrichment (Fig. |4 bottom panels). As observed by |Marm ~ 


Franch et~^ ( |20d9 ), we qualitatively obtain a negative corre- 


lation between galactocentric radius and metallicity, with GCs 
located centrally enriched by up to ^ = 15 polluters, and those 
in the outskirts having ^ = 1-2. 


4. CONCLUSIONS 

The observed age-metallicity distribution of galactic GCs 
shows a bimodal population, with about half the objects re¬ 
siding along an “old” branch with age ^12.8 ±0.6 Gyr and 


spread comparable to the relative ag e error ( [Marm-Franch 
|et al. poo'll Forbes & Bridges|r2010| ). Here, we explored a 


novel idea to form these old GCs at the center of DM miniha¬ 
los with virial temperature T^ir > 10"^ K {M^ ^10^ Mq) and 
with the following additional constraints: 

• Enrichment by metal outfiows from nearby Population- 
II halos rather than by their own Population-Ill star for¬ 
mation; 

• Major merger with a similar metal-enriched, but star- 
free minihalo within a time Af ^ 150 Myr from the 
crossing of the cooling threshold (Ty/r > 10"^ K). 

The minihalo-minihalo merger is the key new ingredient to 
create a compact star cluster. While the details of star for¬ 
mation depend on complex baryonic physics that will be ex¬ 
plored in future investigations, we make robust and falsifi- 
able predictions on the probability distribution of the minih alo 
mergers (age of GCs), which are discussed in Section [3T] We 
obtain {tage) = 13.0 Gyr and /S.tage ^0.2 Gyr, in agreement 
with current observations. 

Improvements in GC age measurements can uncontrover- 
sially falsify the minihalo merger scenario. Also, model pre¬ 
dictions may be within reach of direct observations with the 
James Webb Space Telescope (JWST). In fact, the number 
density of compact star clusters at z > 7 is sufficiently high 
(1.5 X 10"^ arcmin“^) that JWST can identify some of them if 
highly magnified (/i ^ 100) to Mab ^ -17 by a foreground 
galaxy cluster. However, we note that it will be challenging 
to discriminate against othe r compact high-z objec ts such as 
Population-Ill clusters (e.g., |Zackrisson et al.|2012| ). 

Besides GC ages, other aspects of the framework such 
as multiple stellar populations formation triggered by minor 
mergers (Section are more speculative, since they depend 
on a complex interplay between baryonic an d DM physics. 
We presented consistency checks in Section |3.2| to charac¬ 
terize qualitatively the distribution of GCs, their metallicity 
and stripping of the DM halo envelope, but further work is 
required. In addition to detailed investigation of DM strip¬ 
ping, perhaps one interesting aspect to follow-up is the con¬ 
nection between GCs and dwarf galaxies. While earlier stud¬ 
ies considered them as two distinct classes in the luminosity- 
size plane, recent work unveiled the existence of ultra com¬ 
pact dwarf galax ies with intermediate properties (e.g., |Jen-| 
nings et aL]|2014| ). Indeed, our framework predicts a contm- 
uum between classical dwarf galaxies (no major mergers at 
formation, DM dominated) and old GCs (major merger, DM 
envelope stripped by z = 0). Overall, minihalo mergers appear 
to provide a promising scenario to explore for the formation 
of the oldest GCs observed in today’s galaxies. 


It is our pleasure to thank Duncan Forbes, Diederik Krui- 
jssen and Sasha Muratov for useful discussions on an earlier 
version of the manuscript, and an anonymous referee for help¬ 
ful comments. 
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